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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

THEORETICAL PERFORMANCE OF LIQUID HYDRAZINE AND LIQUID FLUORINE 

AS A ROCKET PROPELLANT 
By Sanford Gordon and Vearl N. Huff 


SUMMARY 

Theoretical values of performance parameters for liquid hydrazine 
and liquid fluorine as a rocket propell ant vers calculated on the 
assumption of equilibrium composition during the expansion process for 
four chamber pressures and a vide range of fuel-oxidant and expansion 
ratios- The parameters included were specific impulse, combustion- 
chamber temperature, nozzle-exit temperature, equilibrium composition, 
mean molecular weight, characteristic velocity, coefficient of thrust, 
ratio of nozzle-exit area to throat area, specific heat at constant 
pressure, coefficient of viscosity, and coefficient of thermal conduc- 
tivity. A correlation is presented for the effect of chamber pressure 
on several of the parameters. 

The maximum value of specific impulse was 315.3 pound-seconds per 
pound for a chamber pressure of 300 pounds per square inch absolute 
(20.41 atm) and an exit pressure of 1 atmosphere. 


INTRODUCTION 


Liquid hydrazine and liquid fluorine are of interest as a rocket 
propellant because of high performance. Extensive data exist in the 
literature on their availability, cost, and physical, chemical, and 
handling properties. 


The performance of liquid hydrazine and liquid fluorine has been 
reported in the literature by a number of organizations such as Jet 
Propulsion Lab., C.I.T.; Douglas Aircraft Co., Inc., Project Rand; and 
North American Aviation, Inc. Additional performance calculations for 
this propellant covering a wide range of conditions were made at the 
NACA Lewis laboratory as part of a series of calculations on propellants 
containing the chemical elements hydrogen, fluorine, and nitrogen 
(refs. 1 and 2) to provide a comparison with the performance of other 
propellants based on the same thermodynamic data and computed to the 
same degree of accuracy, and to provide several Ugt^r^ ^ 

viously published. I fill® “ 
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Data vere calculated on the basis of equilibrium composition dur- 
ing expansion for four chamber pressures and a wide range of fuel- 
oxidant and expansion ratios. The performance parameters included are 
specific impulse, combustion-chamber temperature, nozzle-exit tempera- 
ture, equilibrium composition, mean molecular weight, characteristic 
velocity, coefficient of thrust, ratio of nozzle-exit area to throat 
area, specific heat at constant pressure, coefficient of viscosity, 
and coefficient of thermal conductivity. 

A correlation was made which permits the determination of specific 
impulse, characteristic velocity, and ratio of nozzle-exit area to 
throat area for a wide range of chamber pressures. Equations are given 
that permit the calculation of specific impulse for nonisentropic 
expansion or for change in heat content of the combustion gases using 
the originally calculated data. 

Several additional calculations were made assuming frozen compo- 
sition so that data based on the assumptions of equilibrium and frozen 
composition during the expansion process could be compared. 


SYMBOLS 

The following symbols are used in this report: 

A number of equivalent formulas (function of pressure and molecular 
weight; see ref. 3) 

a local velocity of sound, ft/sec 

C p coefficient of thrust 

C° molar specific heat at constant pressure, cal/ (mole) (°K) 
c specific heat at constant pressure, cal/(g)(°K) 


specific heat at constant volume, cal/(g)(°K) 
characteristic velocity, ft/sec 


d log a \ 
3 log T/ s 


( d log pA 
d log T / 
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acceleration due to gravity, 32.174 ft/sec^ 

sum of sensible enthalpy and chemical energy, cal/mole 

sum of sensible enthalpy and chemical energy per unit weight. 


^ n i(0r)i 
J_ , cal/g 


3 


I specific impulse, lb-sec/lb 

k coefficient of thermal conductivity, cal/(sec )(cm)(°K) 

M molecular weight, g/mole 

n number of moles; exponent 

P pressure 

p partial pressure 

q heat, cal/g 

R universal gas constant (consistent units) 

r equivalence ratio, ratio of number of fluorine atoms to hydrogen 
atoms 

S nozzle area, sq ft 

s entropy, cal/(g)(°K) 

T temperature , °K 

w rate of flow, Ib/sec 
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li coefficient of viscosity, g/(cm)(sec) * poise 

p density, g/cu cm 

Subscripts : 

c combustion chamber 

e nozzle exit 

frozen composition assumed frozen 
i product of combustion 

max maximum 

o initial value 

P constant pressure 

s constant entropy 

t nozzle throat 

x any point in nozzle 


0 


CALCULATION OF PERFORMANCE DATA 

Calculations of the performance data were made with a Bell computer 
and an IBM Card -Programmed Electronic Calculator as described in refer- 
ence 1. The assumptions, thermodynamic data, and transport properties 
used for the calculations are the same as those of reference 1. 

The products of combustion were assumed to be ideal gases and 
included the following substances: hydrogen fluoride HF, hydrogen 

H 2 , nitrogen N 2 , fluorine F 2 , atomic fluorine F, atomic hydrogen H, 
and atomic nitrogen N. The dissociation energy of F 2 was taken to 
be 35.6 kilocalories per mole (ref. 4). Physical and thermochemicai 
properties of the propellants were taken from references 3 to 6 and are 
given in table I. 

Procedure for combustion conditions . - The following parameters 
were computed for ten equivalence ratios for a chamber pressure of 
300 pounds per square inch absolute and for five equivalence ratios for 
chamber pressures of 150, 600, and 1200 pounds per square inch absolute: 
combustion temperature, equilibrium composition, enthalpy, mean 
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molecular weight, derivative of the logarithm of pressure with respect 
to the logarithm of density at constant entropy r s , specific heat at 
constant pressure, coefficient of viscosity, coefficient of thermal 
conductivity, and entropy of the combustion products. 

Procedure for exit conditions . - Equilibrium composition, mean 
molecular weight, pressure, derivative of the logarithm of pressure with 
respect to the logarithm of density at constant entropy r s , enthalpy 
of the products of combustion, specific heat at constant pressure, 
coefficient of viscosity, and coefficient of thermal conductivity were 
computed for each equivalence ratio by assuming isentropic expansion 
for four assigned exit temperatures selected to cover the exit pressure 
range from the nozzle-throat pressure to about 0.1 to 0.3 atmosphere. 


Interpolation . - Parameters for pressures at and near the nozzle 
throat, for pressures of 2, 1, 0.5, 0.25, 0.125 atmosphere, and for 
pressures corresponding to altitudes of 10,000, 20,000, 30,000, 40,000, 
and 50,000 feet were interpolated by means of cubic equations between 
each pair of the assigned exit temperatures. The functions and their 
first derivatives used in the interpolations are described in refer- 
ence 1. 


The errors due to interpolation were checked for several cases. 

The values presented for all performance parameters appear to be cor- 
rectly interpolated or in error at most by two or three units in the 
last place tabulated. 

Formulas ♦ - The formulas used in computing the various parameters 
are given in reference 1 and are summarized as follows: 

Specific impulse, lb-sec/lb 

1 ' 294 ' 98 ^To5o^ (1) 

Throat area per unit flow rate, (sq ft)(sec)/lb (pressure in atm) 


s t/ w 


1.3144 T t 

P#t a 


Characteristic velocity, ft/sec 


( 2 ) 


c* * g P c S t /w = 32.174 P c S t /w 


( 3 ) 
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Coefficient of thrust 


C F » Ig/c* » 32.174 l/c* 


Nozzle-exit area per unit flow rate, (sq ft)(sec)/lb (pressure in atm) 


S e / W 


0.040853 T e 
Pe M e I 


Ratio of nozzle-exit area to throat area 

S e/ S t “ (6 

Specific heat at constant pressure, cal/(g)(°K) 

Cp - ^ [t piCCpli +Eni(HT)i Tl *a] 0 

Derivative of the logarithm of pressure with respect to the logarithm 
of density at constant entropy 

5>i D i 

i ft 


Coefficient of viscosity, poise 


FM 

Pi 


Coefficient of thermal conductivity, cal/(sec ) (cm) (°K) 


<=p * ! I) 


When composition is assumed to he frozen, the partial derivatives 
Y. and Y« in equation (7) are equal to zero anu the partial derive- 

tives and in equation (8) are equal to — There- 


fore, equations (7) and (8) become: 


Specific heat at constant pressure assuming frozen composition, 

Cal/(g)( ° K) SntCcJ)! 

(O frozen " ~ Id 
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Derivative of the logarithm of pressure with respect to the logarithm 
of density at constant entropy assuming frozen composition 


(y \ - (°p)frozen ffjA 

'^s 'frozen “ - R/m \c v j fi . ozen 


( 12 ) 


The values of viscosity and thermal conductivity for mixtures of 
combustion gases calculated by means of equations (9) and (10) are 
only approximate. When more reliable transport properties for the var- 
ious products of combustion become available., a more rigorous procedure 
for computing the properties of mixtures may also be justified. 


THEORETICAL PERFORMANCE DATA 

The calculated values of the performance parameters specific 
imnulse, temperature, mean molecular weight, characteristic velocity, 
coefficient of thrust, and ratio of nozzle-exit area to throat area are 
given in tables II to IV. For convenience, the parameters are tabu- 
lated as functions of altitude in table II and as functions of expansion 
ratio in table III. As an aid to engine design, the values of the par- 
ameters within the rocket nozzle for 80, 90, 100, 110, and 120 percent 
of the throat pressure are tabulated in table IV. Equilibrium compo- 
sition, r s , specific heat at constant pressure, coefficient of viscos- 
ity, coefficient of thermal conductivity, and mean molecular weight in 
the combustion chamber and at assigned exit temperatures are given in 
table V. The mole fraction of Fg was not tabulated in table V, since 
it was always less than 0.00003 except at 1600° K for the equivalence 
ratio r = 1.2, where it was 0.00007, 0.00015, and 0.00030 for the 
chamber pressures of 300, 600, and 1200 pounds per square inch absolute, 
respectively. 

Parameters . - The parameters are plotted in figures 1 to 13. 

Curves of specific impulse for various altitudes are shown in figure 1 
plotted against weight percent fuel. The maximum value of specific 
impulse for the sea-level curve in figure l(b) is 315.3 pound-seconds 
per pound at 31.9 percent fuel by weight. 

The maximum values of specific impulse and the weight percentages 
at which they occur were obtained by numerical differentiation of the 
calculated values and sire shown in figure 2 as a function of altitude. 
The maximum specific impulse in figure 2(b) increases 22 percent for 
a chsmge in altitude from sea level to 50,000 feet. 
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Combustion-chamber temperature and nozzle-exit temperature for 
various altitudes axe presented in figure 3 as functions of weight per- 
cent fuel. The maximum combustion temperature in figure 3(b) is 
4411° K at 27.7 percent fuel by weight. The maxi mums of the exit tem- 
perature curves occur near the stoichiometric ratio. 


Characteristic velocity and coefficient of thrust axe plotted in 
figure 4 and the ratio of the area at the nozzle exit to the area at 
the throat is plotted in figure 5 against weight percent fuel. 

Curves of mean molecular weight in the combustion chamber and 
nozzle exit axe plotted against weight percent fuel in figure 6. 

Curves of specific heat at constant pressure, coefficient of vis- 
cosity, and coefficient of thermal conductivity for six pressures axe 
plotted in figures 7 to 9 as functions of weight percent fuel. 

Chamber pressure effect . - The values of the parameters X, c*, 
and S e /S t (or the logarithms of these parameters) are very nearly lin- 
ear with the logarithm of chamber pressure for a fixed equivalence 
ratio and expansion ratio (for example, see fig. 10) . It was therefore 
possible to correlate these data for chamber pressure effects according 
to the following equations to the accuracies indicated: 

I = -^300 (-0*06 percent) (13) 

\" 

J (±0.07 percent) (14) 

S e/ S t = ^ S e/ S t^300 (^) ( ±0 * 3 P ercen *) ( 15 ) 

where I-qq; c* oq , and ( S e/ S t^300 are the values of these parameters 

at 300 pounds per square inch absolute, P c is in pounds per square inch 
absolute , and n is an exponent given as a function of fuel-oxidant 
and expansion ratios for each parameter. The values of n shown in 
figures 11 to 13 were computed from the values of the parameters for 
the various chamber pressures. However, values of n for specific 
impulse can also be calculated from data for a single chamber pressure 
by means of equation (58) derived in the appendix. All the data given 
in this report for X, c*, and S e /Sx are correlated in figures 11 
to 13. 

As an illustration of the use of these curves , suppose it is 
desired to obtain the values of I, c*, and S e /s^ for a combustion 
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pressure of 800 pounds per square inch absolute and an expansion ratio 


of 150 at 30 percent fuel by weight, 
to 13 are: 

I 300 a 380 

°Soo * 7090 

S e /S t - 17.3 

Using these values in equations 
c 800 = 7189 > and (®e/ 8 t^800 “ 16,6 


The values read from figures 11 

n =« 0.0072 
n « 0.0141 
n =* -0.042 

(13) to (15) gives Igoo =* 383, 


Change in heat content of propellant gases . - The results pre- 
sented in this report are computed for adiabatic combustion with pro- 
pellants at the initial temperature indicated. For a change in heat 
content of the combustion gases in the combustion chamber (which would 
result, for example, from heat loss in the combustion-chamber or intro- 
duction of the propellants at a temperature other than the one indi- 
cated), the corresponding specific impulse assuming isentropic expansion 
and the same combustion and exit pressures as in the initial calcula- 
tions may be obtained from the following equation: 

I 2 = I* + B Ah c + C (Ah c ) 2 (16) 

where Ah c is the change in heat content of the combustion gases in 
the combustion chamber. 



and the subscript o indicates the values of the parameters before the 
change is made. The derivation of this equation is given in the 
appendix. 

For example, assume that the enthalpy of the propellant gases fol- 
lowing combustion is 100 calories per gram less than the initial value 
of hj. used in these calculations for the stoichiometric mixture ratio 
at a chamber pressure of 300 pounds per square inch absolute 
(Ahc = - 100). For an expansion ratio of 20.41 (P e = 1 atm), the values 
of I and T obtained from either table III (b) or figures l(b) and 
3(b) and the values of Cp obtained from figure 7(b) are as follows: 
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I 0 = 314.3 lb -sec/lb 
T„ _ = 4394° K 

C jO 

T = 3239° K 

6 jO 

(e ) = 1-958 cal/(g)(°K) 

( c p ) e , 0 = 1 * 122 cal/(g)(°k) 


These values yield the following: 


I 2 = 98,784 
o 

B = 22.87 
C = -0.00278 


By equation 16, 

I 2 = 98,784 + 22.87 (-100) + (-0.00278) (10,000) 

= 98,784 - 2,287 - 28 =.96,469 
I = 310.6 lb-sec/lb 

Equation (16) was applied to a specific case having a Ah c of 
250 calories per gram and resulted in an error of less than 0.1 pound- 
second per pound in specific impulse. For other cases, the relative 
magnitudes of the successive terms in the series may be used to esti- 
mate the approximate magnitude of the error. 


Nonisentropic expansion . - For design of rocket engines and for 
correction of experimental results for heat loss from the nozzle of 
experimental engines, a need exists for performance data assuming non- 
isentropic expansion. In general, the performance assuming nonisen- 
tropic expansion through the nozzle cannot be computed until sufficient 
information is available to determine the heat rejection from the nozzle 
q and the entropy change during expansion As. When suitable informa- 
tion is available” the computation of performance is possible with the 
aid of an enthalpy-entropy diagram or by use of an equation such as the 
following one using data calculated for isentropic expansion: 



87.0132 


+ T q .(As) + 
e * 


This equation is derived in the appendix. 



( 17 ) 
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The use of equation (17) is illustrated by imposing an arbitrary 
relation between q and As by assuming that all the heat q is 
removed from the nozzle throat at constant pressure. For the stoichio- 
metric mixture ratio at a chamber pressure of 300 pounds per square 
inch absolute and an expansion ratio of 20.41, the values of T^. 
obtained from table IV (b) and of (cp)^. obtained from figure 7 lb) are 
as follows: 


T. » 4163° K 

(c p ) t = 1.860 cal/(g)(°K) 

If q » 100 calories per gram, then by equation (44) of the appendix, 


As » -0.02402 - 0.00016 » -0.02418 cal/(g)(°K) 


This value of As, the value assumed for q, and the values of I Q , 

T e ,o> an ^ ( c p)e o 8i ven in the previous example are used in equa- 
tion (17) to give the following result: 

I 2 = 98,784 - 87.0132 jl00 + 3239 (-0.02418) + 2887 ^ .. ♦ . °0 05 q| | 

« 98,784 - 87.0132 (100 - 78.32 + 0.84) = 96,824 

I * 311.2 lb-sec/lb 

Frozen composition . - In order to compare data based on the assump- 
tions of equilibrium and frozen composition during the expansion pro- 
cess, several additional calculations were made assuming frozen com- 
position. These are presented iD table VI together with corresponding 
equilibrium data for the stoichiometric equivalence ratio and for two 
expansion ratios. The percentage differences in these parameters for 
frozen and equilibrium composition are considerably higher for the 
expansion ratio of 163.3 than for the expansion ratio of 20.41. 

For a combustion pressure of 300 pounds per square inch absolute 
and an exit pressure of 1 atmosphere, the values of maximum specific 
impulse are 315.3 pound-seconds per pound at 31.9 percent fuel by weight 
for equilibrium composition during expansion and 292.2 pound-seconds 
per pound at 32.5 percent fuel by weight for frozen composition during 
expansion. 
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Theoretical calculations of the performance parameters of liquid 
hydrazine with liquid fluorine were made for a wide range of chamber 
pressures, expansion ratios, and fuel-oxidant ratios and yielded the 
following results: 

1. For a combustion -chamber pressure of 300 pounds per square inch 
absolute (20.41 atm) and an exit pressure of 1 atmosphere, the maximum 
specific impulse was 315.3 pound-seconds per pound at 31.9 percent fuel 
by weight for equilibrium composition during expansion and 292.2 pound- 
seconds per pound at 32.5 percent fuel by weight for frozen composition 
during expansion. 

2. The maximum combustion temperatures calculated for chamber pres- 
sures of 150, 300, 600, and 1200 pounds per square inch absolute were 
4260°, 4411°, 4539°, and 4687° K, respectively. 

3. The percentage differences in the parameters specific impulse, 
coefficient of thrust, ratio of nozzle-exit area to throat area, exit 
temperature, and exit mean molecular weight due to the assumption of 
frozen or equilibrium composition during expansion are considerably 
higher for the expansion ratio of 163.3 than for the expansion ratio 
of 20.41. 

4. All the data presented in this report for specific impulse, 
characteristic velocity, and ratio of nozzle-exit area to throat area 
have been correlated for chamber -pres sure effect in figures 11, 12, and 
13, respectively. 

5. An equation is given for computing the change in specific 
impulse resulting from a change in the heat content of the combustion 
gases . 

6. An equation is given for computing the value of specific impulse 
for a nonisentropic expansion from data for isentropic expansion. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 9, 1353 



SUMMARY OF RESULTS 
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APPENDIX - DERIVATIONS OF EQUATIONS FOR NONIDEAL PERFORMANCE 
AND FOR EXPONENT IN PRESSURE CORRELATION 

All the data presented in this report were computed on the assump- 
tions of adiabatic combustion and isentropic expansion with propellants 
at the initial temperature indicated. Figure 14 illustrates these 
assumptions on an enthalpy-entropy diagram. In practice, however, 
these ideal assumptions sire not realized. Equations are derived, there- 
fore, which permit using the data in this report to compute specific 
impulse for a change in the enthalpy of the propellant gases (eq. (29) , 
see fig. 15), or for nonisentropic expansion (eq. (37), see fig. 16); 
or to compute the change in entropy for a change in enthalpy along a 
constant pressure line (eq. (44), see fig. 17). 

Correction of specific impulse for a change in enthalpy of propel- 
lant gases. - As may be seen from figure 15, h e is a function of he 
for isentropic expansion from a fixed chamber pressure to a fixed exit 
pressure , 


f(h c ) 


(19) 


Expanding this function in a Taylor's series and neglecting all 
derivatives higher than the second yield 

!>. - f(l>c,o) + f ’< h c,o> + f "< h c,o> ^r- < 20 > 


For constant pressure, 


dh x T da 
dT a JL_ ' 

dh C-n 


( 21 ) 

( 22 ) 


For isentropic expansion. 


ds € 

S„ or "3— 

c ds r 


(23) 


The relations in equations (2l) , (22), and (23) are used to obtain 


f'(hc) 


dh e 

dhl 


TgdSg 

T c ds c 


(24) 


* 
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f"(hc) 


dh,. 


_1_ (^e) 

T c V dh c/ 


+ T e“^r 


2.5^ 

T c dh e dh c 


„ M\ i 1 f!e\ 1 

e \t c 2 / dh c " t c w: w t c 2 cy; 

= !!k_ ["_ — i L_1 (25) 

?c 2 L( c p)e (ep)cj 

Substituting equations (19), (24), and (25) into equation (20) gives 


h = h o + (i? 

e e,o \p 


(26) 

From equation (l). 

I 2 = 87.0132 (h c - hg) 

(27) 

By definition, 

• 



^ “ h c - h c,o 

(28) 

Substituting equations 

(26) and (28) into equation (27) results 

in 


I 2 = 87 


.0132 + bc,o " ^e.o “ ^T c ) 0 

/’ e \ r i i l 

' VT?/o ' W ~J 


By use of equation (27) this may be expressed 


I 2 = I 2 + 87.0132 
o 


L T e\ AV . A 67.0132 , /T c\ f 1 1 1 \ 

l 1 - %i ^ -T- l^il 


t 

( 29 ) 


Equation (29) is identical to equation (16) in the text. 
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Correction of specific impulse for nonisentropic expansion . - As 
may be seen from figure 16, h e is a function of s e for nonisentropic 
expansion, since h g and s g must lie on the line P g . 


f(aj 


(30) 


Expanding this function in a Taylor's series and neglecting all deriv- 
atives higher than the second give 

1*. - *<■..„> + f 'C=e J * ♦ f ^ (31) 


‘e * v e,o 
For constant pressure. 


f(s)=f = T 


f"(s) 


d 2 h dT _T_ 


ds 


,2 “ ds c T 


(32) 

(33) 


Substituting equations (32) and (33) into equation (3l) yields 


»e » V ♦ < T e,o> * + [(^ „ 


(As) : 


(34) 


When equation (27) is modified to account for heat q lost during 
expansion, the result is 

I 2 = 87.0132 (h c - h e - q) (35) 

Taking h c ^ Q = h q and substituting equation (34) into equation (35) 


I 2 - 8^.0132 {h - h^ Q - (T e ^) As 


' _(=p>el 


test . 


(36) 


Equation (27) is used to express this as 


I 2 = I 2 - 87.0132 


f 4 (t ) as + \-^ T \ teg 
V L p C J U 


(37) 


Equation (37) is identical to equation (17) in the text. 

Change of entropy for change o f heat along a constant pressure 
line .“TTs may be seen from figure 17, for a change in enthalpy along 
any line of constant pressure P x , 
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(38) 


Expanding this function in a Taylor’s series and neglecting all deriva- 
tions higher than the second result in 


+ f '(Vo> “ + f " (h x,o> ^ 


(39) 


For constant pressure , 


f’(h) 

dh T 


f’’(h) = L± 


ilk 


-1 dT -1 


dh c 


dh m2 dh 


T c T 


(40) 

(41) 


Equations (38), (40) , and (4l) are substituted into equation (39) to 
give 


s x - s x 


'° + ( T U “ " hi)i 


m. : 


(42) 


When the change in enthalpy at constant pressure P x is due to 
removal of heat q from the system. 


q m -Ah 

(Addition of heat would correspond to negative q.) 

With the aid of equation (43), equation (42) becomes 


(43) 


As = - f — 


(tx/o T^(c ) 

_x'''p'x_|o 


2 


(44) 


An example is given in the text of the use of this equation when 
a nonisentropic expansion is approximated by assuming that the heat is 
removed ax an average pressure. 

Equation (44) may also be used to obtain additional points on the 
constant-pressure lines of an enthalpy-entropy diagram. 


Calculation of chamber pressure correlation exponents . - In the 
general notation, equation (13) may be written 
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Values of n to be used in equation (45) for obtaining specific 
impulse at any chamber pressure P c from the value of specific impulse 
at an initial chamber pressure P c 0 for the same expansion ratio and 
equivalence ratio may be obtained as follows: From equations (27) and 

(45), 

I 2 = 87.0132 (h„-hj = I 2 (Vf (46) 


Differentiating equation (46) with respect to P_ gives 


87 - 0132 (Ife - Sr) = 2n (^) 


Since h c is constant, 


dh * 0 
c 


Substituting equations (46) and (48) into equation (47) yields 


- 87.0132 




87.0132 ( v c 


A /£he\ 
V V dp c/ 


Since entropy s may be expressed as a function of P and h, 




dP c + 


dSe = lap; 


dP e + 


Substituting equation (48) into equation (50) results in 
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ds = ■ 


For isentropic expansion ds e •• ds c j therefore from equations (5l) and 
(52), 


/ds \ /ds \ /ds \ 

14/14/ Pc "U4 


For constant expansion ratio. 


Pc ffc 

p e " dp e 


Substituting equation (54) into equation (53) gives 


9h e\ \( bs o 


e \*TeJ? e Wh c 


dh. ah e \ /a s ; 

dP„ " Iv^pT 


\ 

7h ' 


P e 7 3s c 




Substituting into equation (55) the thermodynamic identity of 
equation (40) and the following identity 


(bs\ = 

WJh 


— - = T 
dP(> e 


k (- ?^)] 


dh~ RT P 


(«c " M e ) 


Equation (57) is substituted into equation (49) to yield 
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The following table presents values of n computed by means of 
equation (58) for the stoichiometric mixture ratio for a chamber pres- 
sure P c o of 300 pounds per square inch absolute. Also presented is 
a comparison of the values of specific impulse at other chamber pres- 
sures computed directly and computed by means of equation (45) from 
these values of n and 1300 * 


hd 

o 

(D 

n 

(from 
eq. (58)) 

P c 
c > 

lb/sq in. abs 

I 

(direct com- 
putation) , 
lb-sec/lb 

I 

(from 

eq- (45)), 
lb-sec/lb 

Error 
in I, 

lb-sec/lb 

20.41 

0.01240 

150 

311.53 

311.58 

0.05 



300 

314.27 


— 



600 

316.92 

316.98 

.06 



1200 

319.45 

319.72 


81.65 

0.00966 

150 

361.27 

361.37 

0.10 



300 

363.80 


— 



600 

366.16 

366.24 

.08 



1200 

368.33 

368.70 

.37 

163.3 

0.00804 

150 

380.48 

380.59 

0.11 



300 

382.72 


— 



600 

384.76 

384.86 

.10 



1200 

386.62 

387.01 

.39 


Equations similar to equation (58) may be derived for evaluating 
the exponent n for c* and S g /S^; however, these equations are 
omitted inasmuch as they involve partial derivatives which have not 
been tabulated in this report. 
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TABLE II. - CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID FLUORINE AS FUNCTIONS OF ALTITUDE 
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^ased on F 2 density of 1.54 at -196° C and N 2 H 4 density of 1.011 at 15 C. 

b Nozzle designed for exit pressure equal to ambient pressure corresponding to altitude indicated. 
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TABLE II. - CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID FLUORINE AS FUNCTIONS OF ALTITUDE - Continued 

(c) Combustion-chamber pressure, 600 lb/sq In. absolute 


4D 


NACA EM E53E12 


25 


I 


Specific 

impulse, 

I, y 

Ib-sec/lb 

LO O O tO tP H 

H CM CM H O <D 
tOTptOcor-r- 

to to to to "O to 

! 1 

*H 0 CD Tp oo to 

<HJ-> 03 CM CO H tO 05 rO 

<H C p Pc, tO LO CO co co r~ 

<D Q) DO I 

O H irHtHHHiHH 

O O O JJ j 

Ratio of 
nozzle- 
exit 
area to 
throat 
area , 

s e /s t 

co 0 CD 

H H to co TP to 
to CM - p r- co to 

10 r- -75 CM r- Tp 
rH rH CM 


CO ^ CO LO OD r- I to -h lo r- co I r- r-- cm cm co 05 


-OlOfOOHH 
tO CO GO 05 O 
tO ro tO tO *0 ^ 


co cm co h cm r- 
to to co r- r- to 


c\j to to co to 

h* to to r* co cn 
to o io 10 o o 


CM LO CO tO r , CM 
fO "B O CO O N 
to to co co r- f- 


LO CO CO to to 
to ^ lo co r- co 
tO tO fO fO to tO 


CT> o r- to to to 

H I s - H CD O ^ 
in lO co co r- r^ 


co cn co r- to OJ 

rH CM tO H* LO CO 

.o to to to to to 


O co o H O 
O ^ o> to r- o 
to co in co co r- 


rH rH rH rH rH rH 


rH LO 

CO ^ 05 O CO 00 
'HT ^ D H H 


CO to OJ 05 05 

cn to m cm >o o co r- r-~ f- r- (M 

o cn to -p h cm r- -p co ro to cn 


tocnr- 

tP CO CM CO CM CO 
tO CO O «H GO tO 


o -C 

c o u, to a> 

(fl rl (j rl E 

0) O H Q) 

s e 3 5 


CO co rH CO CM H 
rH rH CM to 


tO tP CM r- O rH 
I to TP LO LO CO CO 


tO I — O H* O CO 
H H CM CM 


in r- cn to © lo 

rH rH CM 


LO rH LO r- CO CO ] CD GO 05 05 05 05 
O lH rH rH rH rH I H 1 rp M 1 


LO co 05 CM CO to I 
HH CM 


to to *0 to to to 

M 1 M* M* M* M* M 1 


CM CM CM CM CM CM ] CM CM CM CM CJ CM 


OOOOOO GO QO CO CO co GO 
CM CM CM CM CM CM HH rlrlHH 


CD CO CO CO CO CO I 


U - 

05 <V * 
Q. D O 
E 3^0 
<D jo 
Eh aj 


to r- h r- 05 h 
CM rH 05 GO O 
rf CMO r? 

CM CM CM H rH H 


1 Tf CO 05 CO O O 
05 CM tO CM 05 co 
05 03 CO tt H 05 
CM CM CM CM CM rH 


05 to t- go cm 05 I cn o lo r- o 05 ooocMMn 

I -P CO CO tO tP 'O I CM LO CO r" 05 lH 05 Tp 03 H* 05 CO 

CO tP CJ O CO CO [ CM O GO CO tP tO CO LO tO CM O 05 

| CM CM CM CM H <H ! CM CM rH H H rH rH «H rH rH rH 


CO n CM 05 
r- <T> co LO tP 
GO LO 05 GO rH 
O CD TP CM rH rH 


CO tP CO CM 05 
cn co io tt 

CD in <75 XI rH 
O CO TP CM H .H 


o >» 

i H J-> O 
O J-5 *H “ 0) 
<d co o* to 
^ H O 0\ 
«3 D rH X> 

x: <M <D <H 

O -P > 


C Q> D 15) O 

03 rH OJ H E 

Q) O H <y 

E B P 5 


o3 

U 

<D •* «. 

D, <U O ^ 
B U Eh O 
a> p 
Eh X> 


o o o o oo 
o o o o o 
o o o o o 

o o O o o 
H CM -O tP lO 


OOOOOO 

o o o oo 
o o o o o 

* * «k *t 

o o o oo 

H CM to Tp LO 


CO tP CO CM 05 CO Tp CD CM 05 

r- cn co to tp as co lo Tp 

30 tO 05 CD rH i CO LO 05 00 H 

O CO Tp CM H rH 1 O CD tP CM rH H 


co Tp co CM 05 
r- 05 co lO Tp 
CO lO 05 co H 

O CO TP CM H rH 



OOOOOO 
O O O O O 
O O O O O 

o o' a o o' 

rH CM rO Tp to 


OOOOOO 
O O O O O 
O O O O O 
^ ^ ^ ^ 
O O O O o 
H CM to Tp to 


OOOOOO 

O o o o o 
o o o o o 

•> T *k « », 

o o o o o 

rH CM to TP to 


I 

cd * 
> 0) o 

H.OH 
D D 

cr at cd 

W H tr 



ao 

d 


CO 

d 


o 


Based on P 2 density of 1.54 at -196° C and N 2 H 4 density of 1.011 at 15° C. 

‘Nozzle designed for exit pressure equal to ambient pressure corresponding to altitude indicated. 



TABLE II. - CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID FLUORINE AS FUNCTIONS OF ALTITUDE - Concluded 

(d) Combustion-chamber pressure, 1200 lb/sq in. absolute 
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a Based on F 2 density of 1.54 at -196° C and N 2 H 4 density of 1.011 at 15 
^Nozzle designed for exit pressure equal to ambient pressure. 
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a Baaed on F 2 density of 1.54 at -196° C and NgH 4 density of 1.011 at 15° 
^Nozzle designed for exit pressure equal to ambient pressure. 



TABLE III. - CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID FLUORINE AS FUNCTIONS OF EXPANSION 

RATIO - Continued 
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TABLE IV. - CALCULATED PARAMETERS AT PRESSURES NEAR NOZZLE THROAT FOR LIQUID HYDRAZINE 

WITH LIQUID FLUORINE 

{Throat conditions correspond to P x /P t - 1-0. I * (velocity of flov)/gJ 
(a) Combustion- chamber pressure, 150 lb/sq in. absolute 
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TABLE V - PROPERTIES AND COMPOSITION IN COMBUSTION CHAMBER AND FOLLOWING AN ISENTROPIC EXPANSION TO ASSIGNED EXIT TEMPERATURES FOR 

LIQUID HYDRAZINE WITH LIQUID FLUORINE - Continued 

(b) Combustion-chamber pressure, iOJ lb/sq in. absolute « =» 
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(a) Combustion-chamber pressure, 150 pounds per square Inch absolute. 

Figure 1. - Theoretical specific Impulse of liquid hydrazine with liquid 
fluorine. Isentropic expansion to altitude Indicated assuming equilib- 
rium composition. 
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(a) Combustion-chamber pressure, 150 pounds per 
square inch absolute. 

Figure 2. - Maximum theoretical specific impulse and corre- 
spending weight percent fuel in propellant of liquid 
hydrazine with liquid fluorine. Isentropic expansion to 
altitude Indicated assuming equilibrium composition. 
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(c) Combustion -chamber pressure, 600 pounds per 
square inch absolute* 

Figure 2. - Continued. Maximum theoretical specific impulse 
and corresponding weight percent' fuel In propellant of liquid 
hydrazine with liquid fluorine. Xsentropic expansion to 
altitude indicated assuming equilibrium composition. 
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(d) Combustion-chamber pressure, 1200 pounds 
per square inch absolute. 

Figure 2. - Concluded. Maximum theoretical 
specific impulse and corresponding veight 
percent fuel in propellant of liquid 
hydrazine vith liquid fluorine. Isen- 
tropic expansion to altitude indicated 
assuming equilibrium composition. 
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(a) Conibust ion -chamber pressure, 150 pounds per square inch absolute. 

Figure 3. - Theoretical combustion -chamber temperature and nozzle-exit temperature of 
liquid hydrazine with liquid fluorine. Isentropic expansion to altitude indicated 
assuming equilibrium composition- 
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(c) Combustion-chamber pressure, 600 pounds per square inch absolute- 

Figure 3- - Continued. Theoretical combust ion-chamber temperature and nozzle-exit 
temperature of liquid hydrazine with liquid fluorine * Isentropic expansion to 
altitude indicated assuming equilibrium camnositiciu 
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(d) Combustion-chamber pressure, 1200 pound? per square inch absolute. 

Figure 3. - Concluded. Theoretical combustion*^. >' r temperature and nozzle-exit 

temperature of liquid hydrazine with liquid ll*. rine* Isentropic expansion to 
altitude indicated assuming equilibrium composition. 




NACA EM E53E12 



I 




26 3 < 


(a) 

Figure %• - 
hydrazine 
equilibria 
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(b) Combustion-chamber' pressure, 300 pounds per square Inch absolute. 

Fitrure 4. - Continued. Theoretical characteristic velocity and coefficient thrust of 
liquid hydrazine with liquid fluorine. Isentroplc expansion to altitude indicated assum- 
ing equilibrium composition. 
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{ f* Fuel in propellant, percent by weight 

1 ~ 

! (b) Combu3t ion-chamber pressure, 300 pounds per square inch absolute. 

| Figure 5. Continued. Theoretical ratios of nozzle-exit area to throat area for liquid 

| hydrrzlne with liquid fluorine. Isentropic expansion to altitude indicated assuming 

j equilibrium composition. 
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(c) Comb uat ion -chamber pressure, 600 pounds per square inch absolute. 

Figure 5- - Continued. Theoretical ratios of nozzle exit area to th~cs* 

liquid hydrazine with liquid fluorine. Isentropic expansion to altitude indicated 
assuming equilibrium composition. 
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(a) Combustion -chamber pressure, 150 pounds per square inch absolute. 

Figure S - - Theoretical mean molecular weight- in combustion nhamber and at nozzle 
exit for liquid hydrazine with liquid fluorine. Isentropic expansion to altitude 
indicated assuming equilibrium composition. 
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(c) Combustion-chamber pressure, 600 pounds per square inch absolute. 

Figure 6. - Continued. Theoretical mean molecular weight in combustion chamber ard at 
nozzle exit for liquid hydrazine with liquid fluorine. leentropic expansion to alti- 
tude indicated assuming equilibrium composition. 
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(d) Combustion-chamber pressure, 1200 pounds per square inch absolute. 

Figure 6. - Concluded. Theoretical mean molecular weight in combustion chamber and at 
nozzle exit for liquid hydrazine with liquid fluorine. Isentropic expansion to 
altitude indicated assuming equilibrium coinposition. 
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(a) Combustion-chamber pressure, 150 pounds per square inch absolute. 


Figure 7. - Theoretical specific heat at constant pressure of combus- 
tion products (including energy of dissociation) of liquid hydrazine 
with liquid fluorine. Isentropic expansion to pressures indicated 
assuming equilibrium composition. 
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(b) Combustion-chamber pressure, 300 pounds per square inch absolute. 


Figure 7. - Continued. Theoretical specific, heat at constant pressure of combustion 
products (including energy of dissociation) of liquid hydrazine with liquid fluorine. 
Isentropic expansion to pressures indicated assuming equilibrium composition. 
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(c) Combustion -chamber pressure, 600 pounds per square inch absolute. 

Figure 7. - Continued. Theoretical specific heat at constant pressure 
of combustion products (including energy of dissociation) of liquid 
hydrazine with liquid fluorine. Isentropic expansion to pressures 
indicated assuming equilibrium composition. 
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(d) Combust ion- chamber pressure, 1200 pounds per square inch absolute. 

Figure 7. - Concluded. Theoretical specific heat at constant pressure 
of combustion products (including energy of dissociation) of liquid 
hydrazine with liquid fluorine. Isentropic expansion to pressures 
indicated assuming equilibrium composition. 
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(a) Combustion-chamber pressure, 150 pounds per square inch absolute. 

Figure 8. - Theoretical coefficient of viscosity of combustion products 
of liquid hydrazine with liquid fluorine. Isentropic expansion to 
. pressures indicated assuming equilibrium composition. 
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Figure 8. - Continued. Theoretical coefficient of viscosity of combustion products of 
liquid hydrazine with liquid fluorine. Isentropic expansion to pressures indicated 
assuming equilibrium composition. 
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(c) Combust ion -chamber pressure, 600 pounds per square inch absolute 

Figure 8. - Continued, Theoretical coefficient of viscosity of combus 
tion products of liquid hydrazine with licuid fluorine. Isentropic 
expansion to pressures indicated assuming equilibrium composition. 
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(a) Combustion-chamber pressure, 1200 pounds per square inch absolute. 

Figure n. - Concluded. Theoretical coefficient of viscosity of com- 
bustion products of liquid hydrazine with liquid fluorine. Isen- 
tropic expansion to pressures indicated assuming equilibrium 
composition. 
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(a) Combustion-chamber pressure, 150 pounds per square inch absolute. 

Figure 9. - Theoretical coerricient of thermal conductivity of combus- 
tion products of liquid hydrazine with liquid fluorine. Isentropic 
expansion to pressures indicated assuming equilibrium composition. 
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Figure 9. - Continued. Theoretical coefficient of thermal conductivity of combustion 
products of liquid hydrazine with liquid fluorine, Ieentropic expansion to pressures 
indicated assuming equilibrium composition. 
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(c) Combustion-chamber pressure, 600 pounds per square inch absolute. 

Figure 9, - Continued. Theoretical coefficient of thermal conductivity 
of combustion products of liquid hydrazine with liquid fluorine, 
Isentropic expansion to pressures indicated assuming equilibrium 
composition. 
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(d) Combustion-chamber pressure, 1200 pounds per square inch absolute. 

Figure 9. - Concluded. Theoretical coefficient of thermal conductivity 
+ of combustion products of liquid hydrazine with liquid fluorine. 

Isentropic expansion to pressures indicated assuming equilibrium 
composition. 
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Figure 10. - Examp? : iear variation of specific impulse with logarithm of chamber pressure for fixed expan- 

sion ratios. Sto* v .t . „ric equivalence ratio; isentropic expansion to expansion ratio indicated assuming 
equilibrium compos i.v 
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Figure 11. - Theoretical specific impulse for chamber pressure of 
300 pounds per square inch absolute and exponent n for equation 

I = I (P /300) n for liquid hydrazine with liquid fluorine. 
uUU c 

Isentropic expansion to expansion ratio indicated assuming equili- 
brium composition. 
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Figure 12. - Theoretical characteristic velocity for chamber pressure of 
300 pounds per square inch absolute and exponent n for equation 

c* = c^ Q (P^/300) 11 for liquid hydrazine with liquid fluorine. Ise:. 

tropic expansion to nozzle throat assuming equilibrium composition. 
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Figure 13. - Theoretical ratio of nozzle-exit area to throat area for 
chamber pressure of 300 pounds per square Inch absolute and exponent 
n for equation S e /S t - (S e /S t ) 30 0 (Pc/ 300 ) n for lic l uid hydrazine 
with liquid fluorine. Isentropic expansion to expansion ratio indi- 
cated assuming equilibrium composition. 
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